Quantum rings for beginners: 
Energy spectra and persistent currents 
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Theoretical approaches to one-dimensional and quasi-one-dimensional quantum rings with a few 
electrons are reviewed. Discrete Hubbard-type models and continuum models are shown to give 
similar results governed by the special features of the one-dimensionality. The energy spectrum of 
the many-body states can be described by a rotation-vibration spectrum of a 'Wigner molecule' of 
'localized' electrons, combined with the spin-state determined from an effective antiferromagnetic 
Heisenberg Hamiltonian. The persistent current as a function of the magnetic flux through the 
ring shows periodic oscillations arising from the 'rigid rotation' of the electron ring. For polarized 
electrons the periodicity of the oscillations is always the flux quantum <E>o- For nonpolarized electrons 
the periodicity depends on the strength of the effective Heisenberg coupling and changes from <E>o 
first to 4>o/2 and eventually to $o/iV when the ring gets narrower. 

PACS numbers: 73.23.Ra, 73.21.Hb, 75.10.Pq, 71.10.Pm 



C3 



I 

-a 
c 

o 
o 



> 

o 
o 

m 
o 
-i— > 



i 

-a 
c 

o 
o 



I. INTRODUCTION 

The observation of the Aharonov-Bohm oscillations 
and persistent current |2j in small conducting rings, on 
one hand, and the recent experimental developments in 
manufacturing quantum dotsQ and rings 4] with only a 
few electrons, on the other, have made quantum rings 
an ever increasing topic of experimental research and a 
new playground for many-particle theory in quasi-one- 
dimensional systems. 

Many properties of the quantum rings can be ex- 
plained with single-electron theory, which in a strictly 
one-dimensional (ID) system is naturally very simple. 
On the contrary, the many-particle fermion problem in 
ID systems is surprisingly complicated due to enhanced 
importance of the Pauli exclusion principle. In general, 
correlations are always strong, leading to non-fermionic 
quasiparticles as low energy excitations. It is then cus- 
tomary to say that strictly ID systems are not 'Fermi 
liquids' but 'Luttinger liquids' with specific collective ex- 
citations (for reviews see 0, 0, 0, ) • 

The many-particle approach normally used in studying 
the properties of Luttinger liquids starts by assuming an 
infinitely long strictly one-dimensional system. In small 
finite rings a direct diagonalization of the many-body 
Hamiltonian, using numerical techniques and a suitable 
basis, is possible and can provide more direct information 
on the many-particle states. Two different theoretical 
models have been used for the finite rings. In a discrete 
model the ring is assumed to consist of L discrete lattice 
sites (or atoms) with N electrons, which can hop from 
site to site. In this case the electron-electron interaction 
is usually assumed to be effective only when the electrons 
are at the same site. The many-particle Hamiltonian is 
then a Hubbard Hamiltonian 9] or its extension. Another 
possibility is to assume a ring-shaped smooth external 
confinement where interacting electrons move. In this 



continuum model the electron-electron interaction is usu- 
ally the normal Coulomb interaction (e 2 / ^ire^r). In the 
case of a small number of electrons (typically N < 10) the 
many-particle problem is well defined in both models and 
can be solved with numerical diagonalization techniques 
for a desired number of lowest energy states. 

Other many-particle methods have also been applied 
for studying the electronic structure of quantum rings. 
The discrete rings can be solved by using the Bethe 
ansatz^J^J, which becomes powerful especially in the 
case of an infinitely strong contact interaction (the so- 
called i-model). In the case of continuum rings, quantum 
Monte Carlo methods and density functional methods 



have been used (see Sec. 1X111(1 . 

The purpose of this paper is to give an introductory 
review to the many-particle properties of rings with a few 
electrons. Our aim is not to give a comprehensive review 
of all the vast literature published. The main emphasis 
is to clarify the relations between different methods and 
to point out general features of the electronic structures 
of the rings and their origins. Most of the results we 
show in figures are our own computations made for this 
review. However, we want to stress that most of the 
phenomena shown have been published before, (in many 
case by several authors) and we will give reference to 
earlier work. 

We take an approach where we analyse the many-body 
excitation spectrum and its relation to the single par- 
ticle spectrum and electron localization along the ring. 
We will show that, irrespective of the model, the exci- 
tation spectrum in narrow rings can be understood as a 
rotation- vibration spectrum of localized electrons. The 
effect of the magnetic flux penetrating the ring is also 
studied as the change of the spectrum as a function of 
the flux. This is used to analyze the periodicity of the 
persistent current as a function of the flux. Again, it is 
shown that similar results are obtained with the discrete 



2 



and continuum models. 

Throughout this paper we use the term "spinless elec- 
trons" to describe a system of completely polarized elec- 
tron system, i.e. the many-particle state having maxi- 
mum total spin and its z-component. We use lower case 
letters to describe single particle properties and capital 
letters for many-particle properties (e.g. m and M for 
angular momenta). We use terms like "rotational" and 
"vibrational" states quite loosely for describing excita- 
tions, which in certain limiting cases have exactly those 
meanings. When we talk about the rotational state we 
will use the terminology of the nuclear physics and call 
the lowest energy state of a given angular momentum an 
yrast state. 

Experimentally, the study of the spectra of quantum 
ring is still in an early state of an impressive develop- 
ment. It is not yet the time to make detailed comparison 
between the theory and experiments. Nevertheless, we 
will give in Sec. ^a short overview of the experimental 
situation. 

We then attempt to review the theory of quantum rings 
in a logical and pedagogical way, starting with the sim- 
plest case of non-interacting spinless fermions fSecs lIIll 
and llVfl and classical interacting electrons (Sec. [Vj), then 
introducing the effect of magnetic flux (Sec. IVIjl and spin 
fSec lVIl| . Lattice models are presented in Secs lVIlH and 
IIXI followed by numerial approaches (Sees. Ixl IXII IXIIIfl . 
The periodicity properties of the many-body spectrum 
is discussed in Sec lXIII We also briefly discuss the re- 
lation of these previous approaches to the Luttinger liq- 
uid formalism (Sec. IXIV|) . and introduce pair correlation 
functions as a tool to study the internal structure of the 
many-electron state (Sec. IX Vf) . Most of the review will 
deal with rings where the external magnetic flux pene- 
trates the ring in such a way that the magnetic field is 
zero at the perimeter of the ring (Aharonov-Bohm flux) 
and the ring is free from impurities. Nevertheless, in 
Sections IXVTI and IXVTTl we will give short overviews of 
the effects of the Zeeman splitting and impurities on the 
excitation spectrum. 

Several interesting aspects of quantum rings will have 
to be neglected in this paper. Among these are the ex- 
citing possibilities of observing experimentally the spin 
Berry phase (see, however, Sec|n]for a brief experimental 
overview) and, even more exotically, fractional statistics 
(anyons) [T3 |. 



II. EXPERIMENTAL SITUATION 

Since the mid-eighties, there has been an impressive ex- 
perimental development towards smaller and smaller 2D 
quantum rings; with the most recent techniques one has 
reached the true quantum limit of nanoscopic rings con- 
taining only a few electrons pi ll3L IT4T ] . At the same time, 
many of the experiments still study mesoscopic rings with 
hundreds of electrons. Methods of forming such rings in- 
clude litographic methods for forming individual rings on 



a semiconductor surface. The spectroscopic techniques 
are based on the tunneling current through the ring or 
capacitance spectroscopy. Another possiblity is to cre- 
ate a large number of self-organized rings on a substrate. 
The large number of rings allows for observation of direct 
optical absorption. 

One of the hallmarks in this field of research has been 
the experimental observation of the Aharonov-Bohm 
effectjlj or, equivalently, persistent currents. One of the 
main challenges in order to observe this purely quan- 
tum mechanical effect, has been to ensure phase coher- 
ence along the circumference of the ring. Early experi- 
ments in the eighties and nineties reported observations 
of Aharonov-Bohm oscillations a nd p ersistent currents 
in metallic (Au or Cu) rings [TH Il6t and in loops 
in GaAs heterojunctions, i.e. two-dimensional electron 
gas[H, [T^, HD, liH [2^, IlE E3| A related effect which 
has received recent theoretical jH |H ^ ^ and 
experimental|3Cl l3ll l32t l33| attention, is the occurrence 
of a spin Berry phase |34j in conducting mesoscopic rings. 
The simplest example of this topological effect is the 
phase picked up by a spin i which follows adiabatically 
an inhomogeneous magnetic field; the Berry phase is then 
proportional to the solid angle subtended by the mag- 
netic field it goes through. It has been shown [27| that in 
ID rings, a Berry phase may arise due to spin-orbit inter- 
actions. The first experimental evidence of Berry's phase 
in quantum rings was reported in 1999 by Morpurgo et 
al. 30J who interpreted the splitting of certain peaks in 
Fourier spectra of AB oscillations as being due to this 
effect. Very recently, Yang et al.|33j observed beating 
patterns in the Aharonov-Bohm conductance oscillations 
of singly connected rings; these results are interpreted as 
an interference effect due to the spin Berry phase. 

There is by now a vast literature on experiments on 
quantum rings, and in the remainder of this section we 
will just discuss a few selected papers which report mea- 
surements of the many-body spectra, as this is the main 
topic of this review. 

Fuhrer et aL[3(| used an atomic force microscope to 
oxidize a quantum ring structure on the surface of a 
AlGaAs-GaAs heterostructure. By measuring the con- 
ductance through the dot they could resolve the so-called 
addition energy spectrum (see e.g. [33) as a function of 
the magnetic field. The energy levels in the ring and their 
oscillation as a function of the magnetic field could be ex- 
plained with a single particle picture assuming small non- 
spherical disturbation for the ring. The ring had about 
200 electrons. 

Lorke et aL0, 0, [3^ have succeeded to produce a 
few electron quantum rings from self-assembled InAs 
dots on GaAs using suitable heat treating. They 
used far-infrared (FIR) transmission spectroscopy and 
capacitance-voltage spectroscopy to study the ground 
and excited many-body states. The ring radius was es- 
timated to be Rq = 14 nm and the confining potential 
strength Hojq = 12 MeV (the confining potential is as- 
sumed to be \m*io1 l (r — Ro) 2 )- Lorke et al. were able to 



tt 



study the limit of one and two electrons in the ring. The 
experimental findings were in consistence with a single 
electron picture. 

Warburton et al. |39j| studied the photoluminiscence 
from self-assembled InAs quantum rings at zero mag- 
netic field. By using a confocal microscope and taking 
advantage of the fact that each ring had unique charg- 
ing voltages they were able to measure the photoluminis- 
cence from a single quantum ring. The photolumiscence 
spectra show the effect of the Hund's rule of favouring 
parallel spins and other details of the spectra, for exam- 
ple the singlet-triplet splitting. 

The periodicity <f>o = h/e of the persistent current pre- 
dicted for normal metal (not superconducting) rings has 
been observed several times for semiconductor quantum 
rings |l8l Il9l 12(1 l2lL l23j . Also observations of oth er p eri- 
odicities, or higher harmonics, have been reported|24l |40| 
but they have been interpreted as effects of the nonper- 
fectness of the rings. Very recently, the first observation 
of the inherent <&o /N periodicity, which should appear in 
perfect very narrow rings, was reported |4l| . 
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FIG. 1: Configurations of the many-body states of eight 
noninteracting electrons in a ring for spinless electrons. A 
and A' give local minima in the yrast line, while B and C 
correspond to 'vibrational excitations'. 



III. STRICTLY ID-RING WITH 
NONINTERACTING SPINLESS ELECTRONS 

The single particle Hamiltonian of an electron in a 
strictly ID ring depends only on the polar angle <p 



H 



d 2 



2m e R 2 dip 2 ' 



lM*0 = e mv , 



where R is the ring radius, m e the electron mass (or ef- 
fective mass) and mh is the angular momentum. (Note 
that the direction of the angular momentum axis is al- 
ways fixed in a 2D structure). The corresponding energy 
eigenvalues are 



(2) 



2m e R 2 ' 



We will first study spinless electrons, or a polarised 
electron system, where each electron has a S z =T- In 
the noninteracting case the many-body state is a single 
Slater determinant. The total angular momentum is 



M 



and the total energy 



N 



N 



(3) 



(4) 



The lowest energy state for a given angular momentum, 
or the yrast state, is obtained by occupying single parti- 
cle states consecutively (next to each other) . This is due 
to the upwards curvature of e m . We will denote states 



consisting of single particle states with consecutive angu- 
lar momenta, say from mo to mo + N — 1, as "compact 
states" . These states have energy 
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* [Nm lHN*-N)rn + \(2N*- 



(1) and the angular momentum 



M = Nm + 



N(N- 1) 



■TV)] 
(5) 

(6) 



ft is interesting to note that while the single parti- 
cle energy increases with the angular momentum as 
H 2 m 2 /2m e R 2 , the lowest many-body energy increases, 
in the limit of large M, as H 2 M 2 /2(Nm e )R 2 . An yrast 
state of the ring with N electrons corresponds then to a 
single particle with mass Nm e . This is our first notion 
of "rigid rotation" of the quantum state. 

The structure of the yrast states is illustrated schemat- 
ically in Fig. ^ and the actual energy levels as a function 
of the angular momentum are shown in Fig. [5] Naturally, 
the compact states give local minima of the yrast spec- 
trum. In what follows, the most important property of 
the many-particle states is that the internal structure of 
the state does not change when the angular momentum 
is increased by a multiple of N. (By the term "internal 
structure" of a state, we refer to interparticle correlations 
which can be seen by using a rotating frame, as discussed 
by Maksym in the case of quantum dots|42|.') This fol- 
lows from the notion that 

^mTJn({^}) = exp [*!/$>) *S ater ({^}), (7) 



where v is any integer. This kind of change of the total 
angular momentum corresponds to a rigid rotation of the 
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FIG. 2: The lowest total energy of eight noninteracting elec- 
trons as a function of the total angular momentum, i.e. the 
yrast line. Black dots: spinless electrons; open circles: elec- 
trons with spin. The local minima are connected with dashed 
lines. 



state and naturally leads to the above mentioned result 
that the iV-electron system rotates like a single particle 
with mass Nm e . Moreover, correlation funtions are the 
same for both ^f^^jv an d ^M ter since they are derived 
from the square l^l 2 . Note that the minima of the yrast 
spectra occur at angular momenta M — vN if N is odd 
and at angular momenta M = vN + N/2 if N is even. 



IV. LOCALIZATION OF NONINTERACTING 
SPINLESS FERMIONS IN ID 

Classical noninteracting particles do not have phase 
transitions and are always in a gas phase: There is no 
correlation between the particles. In the quantum me- 
chanical case, due to the Pauli exclusion principle, one- 
dimensional spinless fermions behave very differently: 
Two electrons can not be at the same point. This 
means that noninteracting spinless fermions are identical 
with fermions interacting with an infinitely strong delta- 
function interaction. The requirement that the wave 
function has to go to zero at points where electrons meet, 
increases the kinetic energy proportional to 1/d 2 where 
d is the average distance between the electrons (this is 
analogous to the kinetic energy of a single particle in a 
one-dimensional potential box) . The requirement of the 
wave function being zero in the contact points has the 
classical analogy that the electrons can not pass each 
other [3- The pressure of the kinetic energy scales as an 
interparticle energy of the form l/[R 2 (ipi — <fj) 2 }, lead- 
ing to the interesting result that the energy spectrum of 
the particles interacting with the <5-function interaction 
agrees with that of particles with 1/r 2 interaction. In 
fact, the model of ID particles with a 1/r 2 interaction, 



the so-called Calogero- Sutherland modellH El, is ex- 
actly solvable. We will return to this in Sec. I1X El 

We will now demonstrate that the non-compact states 
may in fact be regarded as vibrational states. The sim- 
plest case is that with two electrons. The Slater deter- 
minant is (omitting normalization) 



(8) 



The square of the amplitude of this wave function can be 
written as 



|-0 mi m 2 (<Pi,^2)| 2 = 4 sin 2 



-(Am Ay?) 



(9) 



where Am = mi — mi and Aip = cf2 — This means 
that the maxima of this wave function occur at points 
Aip = (1 + 2n)ix / Am where n is an integer. The com- 
pact ground state has m\ = and mi = 1 occupied, i.e. 
Am = 1 while all noncompact states have Am > 1. This 
means that between and 2n there will be one maximum 
for the ground state, the two electrons being at the op- 
posite side of the ring. For the noncompact states there 
are two or more maxima between and 2ir and the wave 
functions will resemble excited states of a harmonic oscil- 
lator, i.e. vibrational states (actually in the two-electron 
case they are exactly those of a single particle in a ID 
potential box). 

Let us now generalize this analysis to a general case 
with N electrons. To this end we determine the square of 
the many-body wave function, ^({y^})! 2 , in terms of the 
normal modes of classical harmonic vibrations. The equi- 
librium positions of classical particles are ip® — 2nj/N. 
The displaced positions of the particles corresponding to 
a normal mode v are 

^ = ^+^sm(Ki-^27r/A0, (10) 

where v is an integer (1 • • • N/2) and A the amplitude of 
the oscillation. 

Figurc|21shows |^({^ })| 2 for the different yrast states 
of a ring with 8 electrons, as a function of the amplitude 
A of the classical oscillation. In the case of a compact 
state (denoted by A in Fig. ^) |\I/| 2 decreases rapidly 
with increasing A for all v, as seen in the uppermost 
panel. This means that particles appear as being local- 
ized at the sites of classical particles with a repulsive in- 
teraction. For the non-compact yrast states (B and C in 
Fig. njand so on) the maximum of |^| 2 is reached with 
a finite value of A so that v = 1 corresponds to state 
B, v — 2 corresponds to C, etc. The insets of the fig- 
ure show the positions of the electrons at the maximum 
|\I> 2 | of the corresponding figure. Clearly, there seems to 
be a correspondence between the different yrast states of 
noninteracting spinless electrons and vibrational modes 
of classical interacting particles in a ring. 

The conclusion of this section is that the yrast spec- 
trum of noninteracting spinless particles can be under- 
stood as rotational vibrational spectrum of classical par- 
ticles with a repulsive interaction: The compact states 
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FIG. 3: The correlation of different quantum states of nonin- 
teracting electrons with the electron localization at the sites 
determined by the classical harmonic vibrational modes, as 
a function of the amplitude of the atomic displacement from 
the classical equilibrium position, .i.e. ^({v 3 })! 2 as a function 
of A of Eq. ill (H . The panels from top to bottom correspond 

to quantum states (TTTTTTTT), (TTTTTTT o T), (TTTTTT ° TT), 

(TTTTT ° TTT)i (TTTT ° TTTt)i respectively, where o refers to to 
an empty angular momentum state. The solid, long-dashed, 
short-dashed, and dotted lines correspond to v = 1, 2, 3, and 
4, respectively. 
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FIG. 4: Vibrational modes of 40 classical particles in a ID 
ring. Different symbols correspond to different forms of the 
repulsive interaction between the particles, as indicated in 
the figure. The solid line is the result of a nearest neighbour 
harmonic model. The energies are scaled so that the highest 
vibrational energy for each potential is one. 



forward to solve (numerically) the vibrational modes. In 
the case of a short range potential, reaching only to the 
nearest neighbours, one recovers the text-book example 
of acoustic modes of an infinite ID lattice^^ (now only 
discrete wave vectors are allowed due to the finite length 
2ttR) 



are purely rotational states, and the non-compact states 
correspond to vibrational excitations. 



V. CLASSICAL INTERACTING ELECTRONS 
IN A STRICTLY ID RING 

A finite number of classical interacting particles in a 
strictly ID ring will have discrete vibrational frequencies, 
which, after quantization, will give the quantum mechan- 
ical vibrational energies tko v . We assume a monotonic re- 
pulsive pairwise potential energy V(r) between the par- 
ticles, where r is the direct inter-particle distance. The 
potential energy can be written as 



1 N / 



Sill 



(11) 



where R is the radius of the ring and ipi the position (an- 
gle) of particle i. For any pair potential V it is straight- 



fiijj v = C sin 



(£)■ 



(12) 



where C is a contant (proportional to the velocity of 
sound) and v an integer. 

Figure0]shows the classical vibrational energies for dif- 
ferent pair potentials — ln(r), 1/y/r, 1/r, 1/r 2 , and 1/r 7 . 
Naturally, when the range of the potential gets shorter 
the vibrational energies approach that of the nearest 
neighbour interaction, Eci. (|12|) . also shown. Note that the 
vibrational energies of the Coulomb interaction (1/r) do 
not differ markedly from those of l/r 2 -interaction. The 
latter has the special property that the energies agree 
exactly with the quantum mechanical energies of nonin- 
teracting spinless fermions, as discussed in section TlX El 

In addition to the vibrational energy, the classical sys- 
tem can have rotational energy determined by the an- 
gular momentum: E Iot — ^NmR 2 9 (NmR 2 being the 

moment of inertia and 8 the angular velocity). Quantiza- 
tion of the vibrational and rotational energies gives the 
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FIG. 5: Low-energy spectrum of eight electrons interacting 
with infinitely strong 5-function interaction. The black dots 
indicate the state with maximum spin ('spinless electrons'). 
White circles give the energies of other spin states. The spec- 
trum is identical with the rotation- vibration spectrum of eight 
particles interacting with 1/r ^interaction. 



energy spectrum 



VI. EFFECT OF A MAGNETIC FLUX - 
PERSISTENT CURRENTS 

We consider a magnetic flux going through the quan- 
tum ring in such a way that the magnetic held is zero at 
the radius of the ring. This can be modelled, for exam- 
ple, by choosing the vector potential to be (in circular 
cylindrical coordinates) 



^F,if><r c 



2 
B r 



= o 5 -, if r > 

2-rrr 1 



(14) 



which gives a flux $ = irr 2 Bo penetrating the ring in 
such a way that the field is constant inside r c and zero 
outside. If the ring is in the field-free region (R > r c ) the 
electron states depend only on the total flux penetrating 
the ring. 

The solutions of the single particle Schrodinger equa- 
tion of an infinitely narrow ring , 



1 



2m e 



ih d 



55) *»> 



aHlf), (15) 



can be still written as flux-independent plane waves tp m = 
exp(zmy), but the corresponding single particle energies 
now depend on the flux as 



E = E Tot {M) + E yih {{n v }) 



h 2 M 2 

2NmR 2 



(13) 
, the 



where M is the (total) angular momentum and n 
number of phonons v. 

Figure |S] shows the rotational vibrational spectrum de- 
rived from Ea. l|13|) for eight particles with 1/r 2 inter- 
action. It agrees exactly with that calculated quantum 
mechanically for electrons with (an infinite) ^-function 
interaction (see Section 11X0 . It is important to notice 
that since the vibrational states are fairly independent 
of the interparticle interaction, the quantum mechani- 
cal spectrum close to the yrast line is expected to be be 
qualitatively the same irrespective of the interaction. 

Quantum mechanics plays an important role, however, 
when the Pauli exclusion principle is considered. The re- 
quirement of the antisymmetry of the total wave function 
restricts what spin-assignments can be combined with a 
certain rotation- vibration state. For example, for a com- 
pletely polarized ring (maximum spin) only certain rota- 
tional vibrational states are allowed, as shown in Fig. [3] 
Group theory can be used to analyze the possible spin- 
assignments 46, 47] in a similar way as done in rotating 
molecules |48l l49| . 

In the case of bosons (with 5 = 0) the total wave 
function has to be symmetric, and the allowed rotation- 
vibration states are exactly the same as for spin- 1/2- 
particles in the fully polarized state. Figure |3] is then 
also a 'general' result for eight bosons in a quantum ring 
(assuming the interaction to be repulsive). 



e(m, $) = 



2m e R 2 



(16) 



The many-body wave function is still the same Slater 
determinant as without the field, while the total energy 
becomes 



£(M,$) = E(M,0) 



h 2 M 
m e R 2 



h 2 N 



2m e R 2 V $ 



(17) 

The effect of increasing flux is not to change the level 
structure but to tilt the spectrum, say of Fig. |3 such 
that the global minimum of the total energy jumps from 
one compact state to the next compact state. Note that 
Ea. (|17|) is true also for interacting electrons in a strictly 
one-dimensional ring. This follows from the fact that any 
good angular momentum state can be written as a lin- 
ear combination of Slater determinants of non- interacting 
states. For each of these, the last two terms are the same, 
while E(M, 0) depends on the interactions. 

Alternatively, one may perform a unitary transforma- 
tion to obtain a description of the system in terms of a 
field-free Hamiltonian, but with multivalued wavefunc- 
tions ("twisted boundary conditions"). Let us choose a 
gauge A = Vx (which ensures that B = V x A is zero 
at the ring; in a strictly one-dimensional ring one may 
write xiv) — &(p/(2ir) and consider the unitary trans- 
formation 



H - 



H' = UHU' 1 



(18) 
(19) 
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where the operator U is defined as 

JJ — e -ie/Hj Adl _ e ~ie/h x _ 



(20) 



Obviously, the eigenspectrum will be conserved under 
this transformation - if Hip = Eip, then H'ip' = 
UHU^Uip = UHip = Ei//. It is easy to show that 
the effect of this transformation on the Hamiltonian is to 
cancel out the gauge field, i.e. 



U(- 



-iH d „ i — ih d 

- 1 eA 



R dip 



)u- l = 



R dip 



(21) 



and the Hamiltonian takes the field-free form 



H' 



d 2 



2m e R 2 dip 2 



(22) 



Meanwhile, the wave function now picks up a phase 
(—ie/h)x when moving along a given path, even though 
it moves in a region where the magnetic field is zero. This 
is the so-called Aharonov-Bohm effect 0. In particular, 
if an electron moves along a closed path around the flux 
tube, the Aharonov-Bohm phase becomes 



V X • dl = 



2-kR 



A- dl = 



2irR 



-ie 



B • ds = $. 

n 



(23) 



r 



In other words, the boundary condition has now changed. 
While the original wave function satisfies periodic bound- 
ary conditions, ip(<p) — i/j(<p + 27r), for the new 
wave function we have the condition ip'{tp + 2n) — 
ijj'((p)exp(—ie&/h), i.e. "twisted boundary conditions". 
Note that this boundary condition naturally leads to pe- 
riodic eigenvalues e'(m') = e'(m' + $/$o)j m the same 
way as the Bloch condition for electron states in a peri- 
odic lattice leads to the periodicity of the eigenenergies 
in the reciprocal lattice [45|. 

The spectrum l|16|) is clearly periodic in flux with pe- 
riod $0. Any given eigenstate ip' m — exp(i(m— $/<&q)<^), 
however, will have its angular momentum eigenvalue 
shifted by one as the flux is changed by one flux quantum. 

The persistent current of a quantum ring can be writ- 
ten in general as 



/(*) 



dF 
9$' 



(24) 



where F is the free energy of system. To illustrate this 



J 



in the simplest possible case, consider the Schrodinger 
equation for a one-electron ring, 



where 



—D 2 Tp m (ip) = E m ^ m (ip) 
2m 



D _}_d__ie A 

Rdip h 



(25) 



(26) 



Multiplying both sides of Eq.J2HJ by ijj* n from the left 
and integrating along the circumference of the ring one 
gets 



En 



1 



2m 



2tt 



Rd<pip* m (ip)D 2 ip m {ip). (27) 



e JO 



Using the expression l|14|l for the gauge field and taking 
the derivative with respect to flux one obtains, after an 
integration by parts, 



dE m 
d$ 



1 ieh 



2-kR 2m e 



Rdip [ip* m Dip ri 



I 



2tt 



dip j (ip) 



(28) 



where we have identified the RHS with the ID current 
(density) associated with the state m. (Note that j(ip) 
has to be independent of the angle <p). Obviously, the 
same argument applies for a non-interacting many-body 
system where the contributions from different angular 
momentum states are simply summed. The same is true 
in the presence of interactions due to the fact that all flux 
dependence is in the kinetic energy term of the many- 
body Hamiltonian (see discussion after Ea. lfT7|) 1. 



Due to the periodicity of the energy spectrum, the per- 
sistent current will be a periodic function of the flux. In 
the case of noninteracting spinless electrons, the period 
is $0j owing to the fact that the minimum energy for any 
flux corresponds to a compact state. As will be seen in 
Sec. IXIII the electron-electron interactions can change 
the periodicity to <& /2 or to $o/N. 
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When the field is increased, the ground state shifts from 
a compact state (of the kind A in Fig. [fJJ) to the next 
similar state A', i.e. the angular momentum shifts with 
N. At the transition point the three levels, A, D and A', 
are actually degenerate, but the D-states do not change 
the periodicity of the persistent current. However, the 
amplitude of the persistent current is only half of that of 
the spinless electrons. 

So far we have considered only even numbers of elec- 
trons. In the case of an odd number of electrons, the 
spinless case is very similar to that of even numbers, the 
only difference being a phase shift of <&o/2 of the periodic 
oscillations. However, in the case of electrons with spin, 
the odd number of electrons has an effect also on the am- 
plitude of the persistent current: For a large number of 
electrons the amplitude for odd numbers of electrons is 
exactly twice that of even numbers of electrons, as first 
shown by Loss and Goldbart This effect is sometimes 
referred to as a parity effect. 



FIG. 6: Configurations of the many-body states of eight 
noninteracting electrons with spin in a ring. A, A', and D 
give local minima in the yrast line, while B and C correspond 
to 'vibrational excitations'. 



VII. NONINTERACTING PARTICLES WITH 
SPIN 



VIII. NONINTERACTING ELECTRONS IN A 
LATTICE 

Instead of a continuum ring, we will now consider non- 
interacting electrons in a strictly ID ring with a strong 
periodic potential. In this case the standard solid state 
physics approach is the tight-binding model where the 
Hamiltonian matrix can be written as 



The spin degree of freedom allows two electrons for 
each single particle orbital. The yrast states of the many- 
body spectrum are still consisting of compact or nearly 
compact states, but now for each spin component as 
shown in Fig. [f)| The corresponding yrast spectrum for 
8 electrons is shown in Fig. [21 in comparison to the spec- 
trum of spinless electrons. The total spin of the states 
A (in Fig. [fJJ is S — while for all other states have 
either S = (singlet) or S = 1 (triplet) as can easily be 
deducted from Fig. ©by considering the possible ways to 
arrange the ^-components in the orbitals with only one 
electron. 

The yrast spectrum now consists of downward cusps 
at angular momenta M — nN/2, but those minima cor- 
responding to S = compact states are deeper. In the 
case of noninteracting electrons with spin there is not 
such a clear relation to classical rotation-vibration spec- 
trum as in the case of spinless noninteracting case dis- 
cussed above. One could consider spin-up and spin-down 
states separately as spinless systems which do not inter- 
act with each other. However, considering the complete 
many-body spectrum one should notice that the energy 
of state B in Fig. [fj]is degenerate with the state where 
the spin-up system has five electrons and the spin-down 
system only three. 

The persistent current of the noninteracting system 
with the spin is again a periodic function of the flux. 
The period is <i>o like in the case of spinless electrons. 



H, 



eo if j = i 

-t if j = i±l, 



(29) 



where the diagonal terms describe the energy level at 
a lattice site and the off-diagonal terms describe the 
hopping between the levels. This simplest tight-binding 
model assumes one bound state per lattice site and is of- 
ten called the Hiickel model or CNDO-model (complete 
neglect of differential overlap) . 

Assuming the ring, with radius R, to have L lattice 
sites, the problem of solving the eigenvalues becomes a 
text-book example of ID band structure^}, and the sin- 
gle electron eigenvalues can be written as 

e(fc) = e — 2icos(fca) = eo — 2icos ^^~ TO ^ > (30) 

where a is the lattice constant, k the wave vector and m 
an integer. The last step follows from the facts that the 
lattice constant is a — 2ttR/ L, and in a finite ring k will 
have only discrete values k = m/R. Notice that in the 
large L limit, m/L — > 0, the energy spectrum is equiva- 
lent with that of free electrons, Eq. J5J), if we choose the 
tight-binding parameters as eo = 2t and t = h 2 /2m e a 2 . 
In this case the quantum number m gets the meaning of 
the orbital angular momentum. (This equivalence of the 
tight-binding model and free electron model is valid for 
simple lattices of any dimension and can be derived also 
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by discretizing the Laplace operator of the free particle 
Schrodinger equation |5l|). 

The many-body state is again a simple Slater deter- 
minant of the single particle states with the total energy 



E = Ne - 2t cos ( — rrij J , (31) 
j=i V / 

where the selection of the 'angular momenta' rrij is re- 
stricted by the Pauli exclusion principle. For example, 
in the polarized case (spinless electrons) all rrij must be 
different, and the ground state is obtained with a com- 
pact state where the m^'s are consecutive integers as in 
the case of free electrons. Similarly, we can identify 'vi- 
brational states' by making a hole in the compact state 
as will be demonstrated in more detail in the following 
sections. 



113, IH, \^k- It describes electrons on a discrete lattice 
with the freedom to hop between lattice sites, and the 
Coulomb interaction is represented by an on-site repul- 
sion. An interesting feature of the Hubbard ring is that, 
despite being a strongly correlated electron system, it 
can be solved exactly. For a small number of particles 
the solution can be found by direct diagonalization of the 
Hamiltonian matrix as discussed first. For any number 
of particles another solution technique, the Bethe ansatz, 
can be used. This technique is most suitable in the limit 
of infinite U (so-called i-model) and will be addressed 
in the next subsection. The Hamiltonian describing the 
Hubbard model for an iV-electron ring with L lattice 
sites, in the presence of a magnetic flux <j> = <i>/$o pierc- 
ing the ring, can be derived with the help of the unitary 
transformation introduced in Sec I VII ( the flux dependence 
was first derived by Peierls[6(3): 



IX. INTERACTING ELECTRONS ON A 
LATTICE: THE HUBBARD MODEL 

A. Model and exact diagonalization 

A much studied approach to quantum rin gs a nd per- 
sistent currents is the Hubbard model [52L l53l l54l l55l l56l 



N 



H = - t J2H 



N 



n iT n a . 



(32) 



Here, the operator c\ a (ci iCr ) creates (annihilates) an elec- 
tron with spin a at site i; hi a = cj a Ci yCF is the number 
operator for spin-cr electrons at site i. The first part of 
(|32[) describes the hopping of electrons between neigh- 
bouring sites ("kinetic term") while the last part gives 
the repulsion between electrons occupying the same site. 
We will set the hopping parameter t = 1 for simplicity. 

For small numbers of electrons N and lattice sites 
L the Hubbard Hamiltonian can be solved exactly 
by diagonalization of the Hamiltonian matrix. We 
use an occupation number basis (see e.g. H3) 
|*o) = \n a i],n a2 -\ ■ ■ ■n aL f,n a i l ■ ■ -n aL] ) where the z- 
component of the total spin is fixed, i.e. J2i n ai\ = 
and J2i n ail = N l (N = N 1 + N l)- Taking matrix el- 
ements (^oj-ff I'I'q,') gives us a matrix with dimension 
{£ ){]!{ )■ The eigenvalues of this matrix are the exact 
many-body energy levels of the Hubbard Hamiltonian. 
For a given total spin 5", the energy eigenvalues do not 
depend on S z (there is no Zeeman splitting since we as- 
sume the magnetic field to be nonzero only inside the 
ring). In the case of an even N we can choose S z = 0, or 



Nf = Ni = N/2, and the diagonalization of the Hamil- 
tonian will give us all possible eigenvalues (for odd N 
we take S z = 1/2). The easiest way to solve the total 
spin of a given energy state is then to repeat the com- 
putation with all possible values of S z and look at the 
degeneracies. (Note that the matrix dimension is largest 
for S z — and thus solving the matrix for all S z > 
takes less computation than soving the S z — case). 

B. The t- model 

In the following we will first focus on the strong re- 
pulsion limit U — > oo, which was first discussed nearly 
20 years ago|62j: we shall return to finite U effects later. 
In this limit, the system is equivalently described by the 
simpler "t-model" Hamiltonian p3L l64| 

H t = PH km P, (33) 

where Hki n is the kinetic (hopping) term of i|32|l and P 
denotes a projection operator which eliminates all states 
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FIG. 7: Energy spectra of rings of four electrons in five (upper 
panel) and eight (lower panel) sites, calculated with the t- 
model Hamiltonian. The right hand side shows the energy 
levels as a function of the flux (in units of $o). On the left, 
the zero flux energy levels are shown as a function of the 
angular momentum M (black dots). Note that the L — 8 
case shows several vibrational bands. Only negative energy 
levels are shown since the spectrum is symmetric around zero 
energy. 



with doubly occupied sites. It can be shown [63] that this 
projected Hamiltonian is equivalent to a tight-binding 
Hamiltonian describing spinless fermions. Moreover, go- 
ing to large but finite U , near half filling, the Hubbard 
model reduces to the so-called t — J model, i.e. the t- 
model plus a Heisenberg term J £\ Si • |s3 • I n 
other words the translational and spin degrees of freedom 
get decoupled. We will return to the t — J model in Sec. 

EXH 

In the next subsection we will describe how solutions of 
the t-model can be constructed using the Bethe ansatz. 
However, for a small number of electrons and lattice sites 
(N and L) the direct diagonalization has the advantage 
of giving all the eigenvalues at once and some information 
of the many-body state is more transparent. The results 
shown for small discrete rings are obtained either with 
diagonalization techniques or using the Bethe ansatz so- 
lution. It should be stressed that both methods are exact 
and give the same results. 

We will first study four electrons (with spin) in a lat- 
tice having at least one empty site. Figure [7| shows the 
energy spectrum for L = 5 and 8. The right-hand panels 



show the energies as a function of the magnetic flux pierc- 
ing the ring, while the left-hand side shows the lowest 
energies at zero flux, analyzed with respect to the total 
angular momentum, for the different "vibrational bands" 
in the spectra. The following general features should be 
noted: The energy (and thus the persistent current) has 
a periodicity <E>o/W In particular, as first discussed by 
Kusmartsev|54) and by Yu and Fowler 55], in the ground 
state there are N cusps in every flux period; increasing 
the flux makes the ground state jump from one angular 
momentum state to the next. (A more general discussion 
of the periodicity of the ground state will be given in sec- 
tion ^Mj) Each individual state, on the other hand, is a 
harmonic (cosine) function of the flux, with periodicity 
L&q, in accordance with the periodicity of the Hamilto- 
nian. Thus, the energy of each individual state has the 
form, up to an overall phase shift, 



E 



"2tt / 


t $ 


- r -\\ 


T ' 


V $0 


"J. 



(34) 



with p an integer. Obviously, the (kinetic) energy col- 
lapses to zero in the case of half filling, L = N, as in this 
case there is no freedom to hop. Surprisingly, however, 
the above generic flux dependence of the lowest states 
of the spectra is obtained even for a single free site (i.e. 
L = 5 for N = 4). For L > N + 1, there are several "en- 
ergy bands" consisting of states with different amplitude; 
the number of bands increases with the number of empty 

the lowest 



sites. As will be discussed below (Sec. HXDfl 
energy state at any given flux corresponds to the rota- 
tional band without vibrational energy, while the higher 
bands can be interpreted as the vibrational states of the 
system. 

It should be noted that each individual energy level in 
Fig. still has a spin degeneracy: Two or more states 
with different total spin S belong to each energy level. 
The magnetic field can not separate these since we have 
neglected the possible Zeeman splitting. A finite U will 
separate the states belonging to different total spin as will 
be shown in Sec. IIX Fl but will still leave the degeneracy 
due to the S,. 



C. Bethe ansatz 

The Hubbard model in ID can in fact be solved exactly, 
in terms of the Bethe ansatz, as shown by Lieb and Wu in 
1968 11]. The corresponding solution for a Hubbard ring 
in the presence of an Aharonov-Bohm flux has later been 
discussed by a number of authors 0, |§E EE EH- 
Bethe ansatz solutions can be used to construct not only 
the ground state but in fact the entire spectrum presented 
in the previous section. According to this method, the 
energy of a given many-body state may be written as 



N 

cos kj 

3=1 



(35) 



where the numbers kj can be found by solving the set of 
Bcthe equations 
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Lkj = 2nlj - $ - 2tan- 1 [4(sinfc J - - \ f3 )/U] (36) 

0=1 

N JVt 

^2tan- 1 [4(A Q -smkj)/U] = 2nJ a + ^ 2 tan" 1 [2(A Q - \p)/U], (37) 

3 = 1 (3=1 



For obtaining the whole energy spectrum, the unknown 
(complex) constants kj and \p have to be solved for dif- 
ferent quantum numbers Ij and J a , which are restricted 
to be integers or half-integers depending on the numbers 



consider, even N, the Ij:s are integers (half-odd integers) 
and the J a :s are half-odd integers (integers) if iVj is even 
(odd). In practice, all eigenvalues can be found by letting 
p run over all integers and choosing all possible sets of 



N and Nt. For example, if N and Nt both are even, L:s with the restriction I„ 



< L. 



Ij must be an integer and Jj a half-integer. The to- 
tal (canonical) angular momentum corresponding to the 
quantum numbers is given by 



(38) 



a=l 



In the general case of a finite U the nonlinear Bcthe equa- 
tions turn out to be difficult to solve numerically, at least 
for some quantum numbers|67j. For small systems (N 
and L small) it is easier to find all eigenvalues by a direct 
diagonalization of the Hamiltonian matrix as explained 
in the previous section. 

In the limit of infinite U, however, the Bethe ansatz 
solution becomes particularly simple. In this case the 
quantum numbers kj e [0, 27r) are simply given bv |54l l55| 



2?r 



P_ 

N 



$0 



with 



= -2>. 



(39) 



(40) 



Ground state of the t-model 

In analogy to the non-interacting case, the states forming 
the outermost band (see Fig. Q, i.e. states which become 
the ground state at some value of the flux, are compact 
states in the quantum numbers {Ij}. At zero flux, these 
are of course just the yrast states (see left-hand panels 
in Fig. 0. For example, if N and M are both even, the 
quantum numbers {Ij} are consecutive integers, Ij = 
-N/2, -N/2 + 1, N/2 - 1. The corresponding total 
energy, Eq.jnSJ), as function of flux, becomes |55| 



En = -E m cos 



2-7T / $ p 

— — + A 

L I $0 N 



where D r 



(In 



I min) /^ and 



Era 2 



sin(7V7r/L) 
sin(7r/L) 



(41) 



(42) 



So a given solution is constructed for a fixed Nf, i.e. a 
fixed value of the z-component of the spin. The quantum 
numbers Ij and J a are related to the charge- and spin 
degrees of freedom, respectively. In the case we shall at zero flux. Then, Ij 



The integer p should be chosen, for given flux, such as 
to minimize the energy. Of course, there are in general 
several ways of choosing the "spin quantum numbers" J a 
to give the same sum, leading to a large degeneracy of 
the state. This has to do with the fact that spin excita- 
tions are massless in the limit of the i-model; for finite 
U this degeneracy would be lifted. Note that the energy 
collapses to zero for half filling, i.e. when L = N. 

As a simple illustration, consider the example N = 4 
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(43) 

(44) 
(45) 
(46) 



It is easy to check that the "amplitudes", 
— 2 sin(47r/L)/sin(7r/X) agree with those in Fig. 
for L = 5 and 8. Note that the choice of the "central 
value" D c as an overall phase shift makes the sum 
symmetric; one obtains pairs, cos [^5- {—p/N + D c ± 7)] , 
so that, when writing out the sum, all sin(27T7/L)-terms 
cancel. This also implies that shifting all the Ij-.s by 
the same amount, does not alter the total energy 1)410: 
it only changes the overall phase shift D c . The integer 
p is related to the angular momentum M, see Ea. (|38|l . 
This construction gives not only the ground state but, 
for fixed flux, all energy levels of the lowest band, 
corresponding to different p. 
Higher bands 

Generally, excitations can be constructed within the 
Bethe ansatz by introducing holes in the ground state dis- 



tribution of the quantum numbers Ij and J a [6(j . Natu- 
rally, these states are related to the "non-compact" states 
of noninteracting fermions in a strictly ID ring, discussed 
in Sec. IIIII The higher bands in our spectra correspond 
to charge excitations, i.e. holes in the Ij ( J Q -excitations 
do not cost any energy in the infinite U limit), and this 
reproduces exactly all the energies obtained numerically 
in the previous section. We will illustrate this procedure 
with a few examples for N = 4, M = 2. The lowest 
possible "excitation" in the charge quantum numbers Ij 
is lifting the topmost one by one step, 

{Ij} = -2, -1,0, 2, (47) 

and the corresponding kj as given in Ea. H39|) . The re- 
sulting energy is 



Ei = — 2^ cos kj 
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(48) 



(49) 



This again gives a band of states which are cosine functions of the flux, with period L$o- Differentiating Ei wrt. 
a = ^ (-■£■ - - 

u — L \ N 



one finds the "amplitude" of this band as 
Ei, 



1 n 47r 27r 1 

mm = —2 1 + 2 COS — + COS — COS Ot mln 



L 



L 



n ■ 27r • 

2 sm — sm a mm 



(50) 



with 



= tan 



sin^ 



_ 1 + 2 cos ^ + cos ^ 



2tt 



(51) 



The simplest example for four particles is L = 6. In 
this case a m i n — 7r/3, i.e. the energy of the state 1471) has 
a minimum at zero flux with p/N = — 1 and -Bi, mm = 
— cos(7r/3) — v3sin(7r/3) = —2. 



Note that at zero flux, the lowest state of this band may 
have an energy larger than Ei >m i n - the nearest minimum 
may occur at a finite flux. This is because p has to be 
an integer, and (LN/2ir)a m i n is not generally an integer 
(except for some special values of L) . 



This method is easily generalized to construct higher 
bands. For example, the second excitation band corre- 



sponds to, for four particles, {Ij} 



-1,1,2 with 
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TABLE I: Total weight of the most important configurations 
of the many-body states of different vibrational states for the 
t-model ring with L = 8 and N = 4. The configuration 
is shown as filled and empty circles indicating whether or 
not there is a electron in the corresponding site. The second 
column shows the number of such states (not including the 
degeneracy coming from the spin configuration). The last 
three colums show the total weights of these states for the 
vibrational ground state (wq) and the two vibrational bands. 
The total weights are the same for all rotational states and 
different spin states belonging to the same vibrational mode, 
configuration n wo wi u>2 
o • o • o • o» 8 0.1248 0.0000 0.0000 
o • o • • o »o 32 0.1824 0.2136 0.0624 
o • • o o • mo 16 0.0624 0.0000 0.1256 
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.(52) 



The higher (inner) bands constructed in this way may 
be interpreted as corresponding to vibrational excita- 
tions. As an example, we have examined, for the N = 4, 
L = 8 solutions of the t-model, which electron-hole con- 
figurations have the largest amplitude (in the many-body 
wavefunction) in each of the three lowest bands. Table 
[J] shows the weights of the most important basis states 
for the ground state and for the two lowest vibrational 
states. These are consistent with the classical motion of 
electrons in the first vibrational modes. 



given by E n 



J2iLi m i ~ M 2 /N. As discussed in Sec. 



11111 the ground state corresponds to filling consecutive 
angular momentum states, e.g. rrii — 0, 1, ...,N — 1. A 
set of N/2 excited states is constructed by creating 1-hole 
excitations, the lowest one being rrii — 0, 1, N — 2, N. 
(Note that an overall shift of all the angular momenta 
does not change the minimum energy - it just occurs at 
a different flux.) Computing the ratios R n as defined in 
Eg. (15*3*1) one finds 



Ri 



_ n(N - n) 
= N - 1 



(55) 



As we shall see, these one-hole excitations in a sense cor- 
respond to fundamental phonon excitations. Moreover, 
it is easy to show that 2-hole excitations lead to ratios 
(in the limit N — » 00) which are twice those in Eq.JSSJ, 
thus corresponding to two-phonon excitations, etc. 

Now, the same ratios can be recovered from the Bethe 
ansatz solution of the t-model, as given by equations (|35|> 
and (|39|l : Constructing the set of excitations which corre- 
spond to one hole in the quantum numbers Ij , minimizing 
the energy as described in the previous section, and tak- 
ing the limit L — > 00, one again gets Ea. (|55|l . This shows 
that the energy bands of the Hubbard model reduce to 
those of non-interacting particles in the limit of infinite 
repulsion and infinite number of sites. (Of course the cor- 
respondence between the t-model and free particles can 
also be seen using a similar argument as in section "Villi 
For the maximum spin state, — p = J a = in Ea. l|38() 
so that the integers Ij correspond just to the single par- 
ticle angular momenta.) 



D. Vibrational bands 

In the limit L — > 00, i.e. infinitely many sites, the 
t-model is expected to correspond to a system of non- 
interacting (spinless) particles in the continuum (as they 
no longer "see" the delta function interaction) , in a sim- 
ilar way as the simple tight-binding model approaches to 
the continuum model when the number of lattice sites 
increases (see Sec. IVIIII) . One way of illustrating that 
this is indeed the case, is to examine the ratios 



E n — Eq 

Ei — E 



(53) 



where Ei is the (minimum) energy of the ith excited 
band. In the case of non-interacting, spinless particles 
the energy, in units of h 2 / (2m e R 2 ) , of an A^-particle state 
at flux cj) = $/"I>o is given by 



E = J2(r 



(54) 



where rrii are the single particle angular momenta and 
the total angular momentum is M — ^2 rrii. The mini- 
mum in energy of this state occurs at = M/N and is 



E. 1/djj interactions: 



There is another type of interaction in one dimension 
which leads to the same set of excitation bands as above, 
namely Vij = l/df^, where dij is the distance between 
particles i and j. Models in ID with this type of inter- 
action are known as the Calogero- (on a line) or Suther- 
land model (on a circle) [43L and have been studied 
extensively over the past three decades. What makes the 
Calogero (or Sutherland) model special, is that it mimics 
as closely as possible a system of free particles. Consider 
the A^-particle Sutherland Hamiltonian, 



N 



^ = -Ea? + 2A ( A - 1 )E 



^(2sin(^))- 



,(56) 



where we have set the radius of the circle to 1. The 
complete excitation spectrum of this model can be found 
exactly in terms of the asymptotic Bethe ansatz (see e.g. 
|68| and references therein) and is remarkably simple. 
The total energy and angular momentum of a given state 
can be written as E = J^jP'j an d P = J2jPji respec- 
tively, where the quasi-momenta p n are related to the 



14 



free angular momenta rrij by 



Pi = m j + 7\J 



N+l 



(57) 



where 7 = (A — 4)/4. In other words, the energy and 
total angular momentum take the same form as in the 
non-interacing case, with the interactions absorbed in 
the shift of the quasi-momenta. For each non-interacting 
many-body state characterized by a set of fermionic 
quantum numbers {rrij} there is a corresponding state 
with quantum numbers {pj}. Note that P = ^2jPj = 
J2j m j = M. Now consider the Sutherland model in the 
presence of an Aharonov-Bohm flux piercing the ring. 
The energy of a given state can then be written as 



(58) 



and, as in the non-interacting case, is minimized for cf> — 

M/N, 



3 



p) - M 2 /N. 



(59) 




We want to compute the shift in energy of a many-body 
state as compared to the non-interacting one, which then 
will be used to compute the ratios R n . Inserting the solu- 
tion H57|) into Eg. (|58|l one finds, up to an overall constant, 

AE{j) = £ mm ( 7 )-£ mm ( 7 = 0) (60) 



jrrij - (N + l)M . (61) 



Again constructing the set of one-hole excitations (in rrij) 
as before, one easily finds the difference between the shift 
of the nth excitation and the shift in the ground state, 

A£W (7) - AE^ (7) = 7 (n{N - nj) . (62) 

From this it immediately follows that the ratios R n 
for the Sutherland model are the same as in the non- 
interacting case, see Eq. l|55|l . In addition, recall section 
[V]where we performed a calculation of the normal modes 
for a set of classical particles on a ring with 1/d 2 repul- 
sion. In a semiclassical picture, the corresponding fre- 
quencies ujj were then interpreted as the eigenfrequen- 
cies of the many-body problem, with energies fvjjj. Us- 
ing these to compute the ratios R n , one again obtains 
the same expression as above. This illustrates that the 
"single-particle excitations", both in the non- interacting 
case and in the Sutherland model, can indeed be inter- 
preted as corresponding to vibrational modes. (Cfr. the 
discussion of the free particle case in Sec. II VI ) 



F. Finite U 

The Hubbard model with finite U can still be solved 
with the Bethe ansatz but now it requires numerical so- 
lution of the set of nonlinear equations i|36|) ■ In the case 



N=4 
L=4 



FIG. 8: Energy levels as a function of U for a four electron 
Hubbard ring with four sites. 



of small number of electrons and sites a direct numeri- 
cal diagononalization of the Hubbard Hamiltonian, Eq. 
(|32[1 . is in fact easier. The results shown in this section 
have been computed with the direct numerical diagonal- 
ization. As an example case we use again the 4 electron 
ring. 

Figure [S] shows the energy spectrum for four electrons 
in four sites as a function of on-site interaction U . The 
results start from the noninteracting case (with spin) . In- 
creasing U separates the spectrum into different groups. 
The lowest group corresponds to states where all the elec- 
trons are mainly at different lattice sites while the two 
higher bands correspond to states where one or two elec- 
trons are at the same site, respectively. This can be seen 
by looking at the structure of the many-body states or 
simply by noticing that the energy of the higher groups 
of levels increase as E as U and E as 2U. 

The high-C/ limit of the lowest group of states can be 
explained with the so-called t J-model. In the limit of 
large U the half-filled Hubbard model can be approxi- 
mated as (see [H3 an d references therein) 



H 



T N 

-Y 

2 ^ 



Sj ' S-j 



(63) 



where H t denotes the Hamiltonian of the infinite U limit 
(t model), J = 4t 2 /U and Si is a spin operator (S = 1/2). 
In the case of a half-filled band the t-model gives only 
one state (with zero energy). In the large U limit all 
the low energy states can then be described with an an- 
tiferromagnetic Heisenberg Hamiltonian which separates 
the different spin-states. The case with four electrons is 
especially easy to solve (see exercise 30.3. of Ref. |45(), 
leading to energy levels shown in Table[H] The table also 
shows the orbital angular momentum determined from 
the symmetry properties of the Heisenberg state. 

In the case of a large U the magnetic flux does not 
have any effect if N = L due to the fact that the electron 
motion is strongly hindered. The situation changes if 
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TABLE II: Energy states of a Heisenberg ring with four elec- 
trons, sorted according to the spin S and total angular mo- 
mentum M quantum numbers. M is determined with the 
help of the symmetry of the state. 

S S z ME 






-1,0,1 
-1,0,1 
-1,0,1 
-1,0,1,2 



-2J 
-.1 




J 




FIG. 9: Flux and U dependence of the many-body states of 
the Hubbard model with N = 4 and L — 8. For each fixed U 
the flux goes from to $0; in between, $ = and U increases 
linearly to the next fixed value. The lowest overall energy and 
the lowest energy state corresponding to the maximum spin 
are shown as thick lines. Note that the maximum spin state 
is independent of U and that the periodicity of the yrast state 
changes from $0 to $o/4 when U increases from to 00. 



empty sites are added. Figure shows the development 
of the low energy states as a function of the magnetic field 
and U in the case of four electrons in eight sites. When U 
is reduced from infinity, the different spin-states separate, 
causing the periodicity of the yrast line to change from 
^o/TV at U = 00 to $ at U = 0. This cha nge o f the 
periodicity is addressed in more detail in Sec. IXIII 

The effect of finite, but large U, is to split the degen- 
eracy of the different spin-states of the i-model. If there 
are no empty sites, the Hamiltonian approaches that of 
the iJ-model, Eq. (|u3|) . We have solved the spectrum 
of four electrons as a function of the flux, increasing the 
number of sites from 4 to 12. The results show that when 
the number of empty sites increases, the spectra are still 
in fair agreement with those of the t J-Hamiltonian, but 
the effective coupling between the spins, J, is not any 
more At 2 /U , but decreases rapidly when the number of 
empty sites increases. 

Figure ^| shows the difference of the energy spectra 
derived from the exact Hubbard spectrum and from the 
t-model, for different values of empty sites. The results 
are scaled so that the energy difference between the two 
spin states corresponding to M = are the same. We can 



-1 



-3 



12 3 

angular momentum M 

FIG. 10: The Heisenberg model energy spectrum for four 
electrons determined from the Hubbard model with empty 
sites. The black dots are the results for four electrons with 
U — 100 and different number of sites: For each state the 
dots from left to right correspond to L = 5, 6, 8, 10, and 
12, respectively. For comparison, the solid lines denote the 
half-filling case (no empty sites), i.e. L — 4. In each case the 
energy difference of the two M = states is adjusted to be 1. 
Open circles show the same for U — 10. 



2 4 6 8 10 12 14 

L 

FIG. 11: The dependence of the effective Heisenberg coupling 
J on the number of sites (L). The crosses show the asymptotic 
large U limit of Eq. Iti4ft . the black dots (open circles) the 
result determined from the exact solution for U = 100 (U = 
10). 



see that when the number of empty sites increases, the 
difference spectrum approaches that of the Heisenberg 
Hamiltonian. These results suggest that the separation 
of the Hubbard Hamiltonian (for large U) into the t- 
model and Heisenberg Hamiltonians is accurate in the 
two limits, N/L — > 1 and N/L — > 00. Moreover, the 
agreement with all numbers of empty sites is surprisingly 
good, even if U is as small as 10. 

Yu and Fowler |55j have used the Bethe ansatz to study 
the large U limit for any number of empty sites and shown 
that the effective J c ff is related to the quantum numbers 
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3=1 



Figure ^] shows the effective coupling constant deter- 
mined by this equation as a function of the number of 
sites for a four electron ring, compared to those obtained 
from the direct diagonalization of the Hubbard Hamilto- 
nian for U = 10 and U = 100. The agreement is fairly 
good for any number of empty sites even for the relatively 
small U — 10. The decrease of J e g as a function of L can 
be explained with electron 'localization': When the num- 
ber of empty sites increases, the localized electrons move 
further apart reducing the exchange interaction. We find, 
once again, that in the ID system the electrons with re- 
pulsive (even ^-function) interactions behave as localized 
particles. 

X. QUASI-1D-CONTINUUM RINGS: EXACT CI 
METHOD AND EFFECTIVE HAMILTONIAN 

In this section we review electronic structure calcu- 
lations for quasi-lD rings, published in Refs. |46l l47| : 
they were a natural continuation of earlier calculations 
done for electrons in harmonic two-dimensional quantum 
dots|37|. Usually the quantum ring is described with a 
displaced harmonic confinement (although several other 
models have been used [63, rTPj l 

V(r) = ^m e u; 2 (r-r a ) 2 , (65) 

where r is the radius of the ring and loq the perpendic- 
ular frequency of the ID wire. Note that we still assume 
the ring to be strictly two-dimensional, i.e. it is infinitely 
thin in the direction perpendicular to the plane of the 
ring. The parameters describing the ring, ro and loq can 
be related to the density parameter r s — l/(2no) of the 
ID system (no is the ID density) and to a parameter de- 
scribing the degree of one-dimensionality of the wire. For 
the latter, Reimann et al. |46l FH used a parameter C f 
defined with the relation fkj = Cf^tt 2 / '(32mr 2 ). The 
physical meaning of Cf is that it is the ratio of the first 
radial excitation in the ring to the Fermi energy (approx- 
imated by that of an ideal ID Fermi gas with the same 
density) of the ID electron gas. 

There are several approaches to solve the many-body 
problem of interacting electrons in the above poten- 
tial. In studying the spectral properties the most useful 
method is the brute force diagonalization of the many- 
body Hamiltonian in a proper basis. This Configuration- 
Interaction (CI) method gives the whole many-body en- 
ergy spectrum as well as the corresponding many-body 
states. Naturally, the solution can be numerically accu- 
rate only for small numbers of electrons, typically less 
than 10. The matrix dimension can be reduced by fix- 
ing the orbital angular momentum to the desired value. 




M=0 Angular Momentum M *■ M=6 
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FIG. 12: Energy spectra for two quasi-one-dimensional con- 
tinuum rings with six electrons (in zero magnetic field). The 
upper panel is for a narrow ring and it shows several vibra- 
tional bands. The lower panel is for a wider ring which shows 
stronger separation of energy levels corresponding to differ- 
ent spin states (shown as numbers next to the energy levels). 
Note that also the narrow ring has the same spin-ordering of 
the nearly degenerate state as expanded for the lowest M = 
state. 

For example, Koskinen et al. 46] first expanded the solu- 
tions of the single particle part of the Hamiltonian in a 
harmonic oscillator basis and then used these functions 
as a single particle basis for the Fock space in doing the 
CI calculations. According to their eigenenergies, up to 
50 single-particle states were selected to span the Fock 
space and the number of the many-body Fock-states was 
restricted to about 10 5 using another energy cutoff (for a 
given total M). For an even number of electrons all spin- 
states can be obtained with fixing Nf = N± = N/2. The 
total spin of each state can afterwards be determined by 
calculating the expectation value of the S 2 operator. 

FigurclT^lshows the calculated energy spectra obtained 
with such a calculation for two different rings with six 
electrons in each. It is instructive to introduce a simple 
model Hamiltonian |46j 

h 2 

Hef f = jjM 2 + + J Si • Sj (66) 

where uj v are the frequencies and the integers n v the num- 
ber of excitation quanta of the various vibrational normal 
modes, and I is the moment of inertia of the 'molecule'. 
This Hamiltonian is thus simply a combination of rigid 
rotation of the whole system, internal vibration, and a 
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FIG. 13: Low energy spectra of a four electron quantum dot 
(upper panel) and a four electron quantum ring (lower panel) 
as a function of the angular momentum. The numbers next 
to the energy levels give the total spin. CI calculation from 

Heisenberg term to capture the spin dynamics, and one 
may examine how well it describes the exact results. To 
this end, note the following interesting features in Fig. 
1121 The narrower ring shows clear rotation-vibration 
bands, very similar to those obtained for electrons in a 
continuum ring with 5-function interaction, cfr. Fig. [S] 
The only difference is that now the ratios of the vibra- 
tional states correspond to those determined for classical 
electrons interacting with 1 jr interaction. 

Each spectral line consists of several nearly degener- 
ate spin-states. The inset shows as an example the de- 
tailed structure of the M = state. The spin structure 
coincides with that determined from the t J-model, i.e. 
it can be determined by solving the antifcrromagnetic 
Heisenberg model for a ring of six electrons. The ratios 
of the energies of different spin-states are quantitatively 
the same as in the Heisenberg model. 

The lower panel of Figure ^] shows that this is true 
also for a wider ring. In this case only the vibrational 
ground state is clearly separated from the rest of the 
spectrum. However, the internal structure of this yrast 
band is still very close to that of the Heisenberg model: 
Qualitatively the agreement is exact, i.e. each angular 
momentum has the right spin states in the right order. 
Only the energy ratios are not any more exactly the same 
as in the Heisenberg model. Koskinen et al.|4j] have 
studied in more detail how well the model Hamiltonian 
(|66fl describes the exact many-body results for a ring of 
six electrons. 



Exact CI calculations for a ring of 4[72|, and 5 and 
7|4^ electrons show similar agreement. Figure IT31 shows 
the energy spectra for a four electron ring with a compar- 
ison to a four electron dot. With such a small number 
of electrons even the dot shows nearly the same yrast 
spectrum as the ring. This indicates that also in the dot 
the electrons will localize in a square Wigner molecule 
(see also jzjj). Similar localization of electrons and their 
rotational and vibrational spectra are observed also in 
two-electron quantum dots l74l (where the energy spec- 
trum can be solved exactly 75]). When the number of 
electrons in a quantum dot is 6 or larger the classical 
configuration of electrons in the Wigner molecule[7(| is 
not any more a single ring and the spectral properties 
become more complicated. Nevertheless, even there the 
polarized case shows rotational bands consistent with the 
localization of electrons jx^]. We should point out that the 
idea of describing few-electron systems in terms of rigid 
rotation and internal vibrations, is not new, but was first 
applied to quantum dots by Maksvm[42^. 



XI. EXACT DIAGONALIZATION: FINITE 
MAGNETIC FIELD 

Flux inside the ring 

Let us first consider a magnetic field concentrating in 
the center of the ring, as described with the vector poten- 
tial (d|. Since the r-component of the vector potential 
is zero and the 93-component is proportional to 1/r, the 
only effect of the flux is to change the angular momen- 
tum term of the single particle Schrodinger equation as 

ft 2 m 2 ^ fi 2 (m-$/$ ) 2 
2m e r 2 2m e r 2 

It is easy to see that for integer $/<!>o the energy levels 
and single particle states are equivalent to those without 
the flux, they are just shifted to another angular momen- 
tum value. Since each Slater determinant of the many- 
body state has a good angular momentum quantum num- 
ber, the same is true for the many-body state: The angu- 
lar momentum M is shifted to that of M — Nv when a in- 
teger number v = <&/&q of flux quanta is penetrating the 
ring (as in Eq. (|17fl for a strictly ID ring). This means 
that we get the same picture as discussed before: The 
effect of the increasing flux is just to tilt the spectrum so 
that successively higher and higher angular momentum 
values become the ground state. The ground state energy 
will be periodic with flux: E M {§) = E M +Nv{$ + v&o)- 
Note that this is true even if the M-dependence of the 
many-body spectrum is not exactly proportional to M 2 ; 
the only requirement is that the flux is restricted in the 
central region of the ring and the magnetic field does 
not overlap with the single particle states. In this case 
we can write the model Hamiltonian describing the yrast 
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spectrum of the exact CI calculation as 

h = j • Sj 



2Nm e R 2 



M- — \ . (68) 



The observed periodicity of the ground state energy, or 
persistent current, as a function of flux is determined by 
the variation of the yrast-line as a function of M, for 
example as shown in Fig. El If the situation is like 
in the upper panel of Fig. 1121 i.e. when the ring is 
very narrow, the energy increases accurately as M . In 
this case Eq.(|B^|l gives the periodicity $o/-W- On the 
other hand in the case of the lower panel of Fig. IT21 
the minimum energy jumps from M = to M = 3 and 
then to M = 6, when the flux is increased. This means 
a periodicity of $0/2- If the ring is made even wider, 
eventually the minimum at M — 3 will not be reached 
and the periodicity changes to <f>o- 

Finally, we should notice that if the electron gas were 
polarized (spinless electrons), the periodicity would al- 
ways be $o- This can be seen from Fig. 1121 which shows 
that the maximum spin state occurs in the lowest vibra- 
tional band only at the angular momentum M — 3 (more 
generally, at M = N/2 for even number of particles and 
at M = for odd number of particles) . 



Homogeneous magnetic field, no Zeeman splitting 

Experimentally it would be easier to measure the quan- 
tum rings in the presence of a homogeneous magnetic 
field. This case was first treated by first principle calcu- 
lation methods in the pioneering papers by Chakraborty 
and PietilainenjT^ and by Niemela et al|7(i|. and the 
present section is mainly based on these papers. In this 
case the vector potential can be expressed, for example, 
in terms of a symmetric gauge A = ^(—By, Bx, 0). This 
vector potential effectively adds an additional harmonic 
confinement centered at the origin. The r-dependent sin- 
gle particle potential changes as 



1 



m e ul(r -r y 



1 



m e uj B (r — re) 2 + constant (69) 



where uj 2 b = loq + e 2 B 2 /{Am 2 ) and r B — t uj 2 /uj 2 b . The 
effect of the field is then just to change the parameters 
of the confining ring. This will change the energy differ- 
ences between the single particle states, but as long as 
the ring is narrow enough to have only one radial mode, 
it can not change their order. Consequently, in narrow 
rings the effect of the field on the many-body state is 
small and only quantitative. Nevertheless, the change of 
the potential shape and the constant term means that the 
lowest single particle state increases with the fluxjT^, as 
shown in Fig. Elf° r spinless electrons. The figure also 
indicates that the effect of the electron-electron interac- 
tions in the spinless case is mainly to shift the spectrum 
upwards by a constant. 
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FIG. 14: Energy levels for four noninteracting (dashed lines) 
and interacting (solid lines) spinless electrons in a ring as a 
function of the magnetic flux. The magnetic field is homoge- 
neous. From Ref. |7q|. 
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FIG. 15: Energy levels as a function of the magnetic flux in 
a ring of four electrons. The left-hand side shows the results 
for delta-function interaction in a strictly ID ring, (a) for 
noninteracting, (b) and (c) for interacting so that in (c) the 
interaction is twice as strong as in (b). The right hand panel 
shows the result for a quasi-lD ring with 1/r- interaction, (c) 
noninteracting and (d) interacting electrons. Note the simi- 
larity of the results in (b) on the left-hand side and (d) on the 
right hand side. From Ref. |7fJ 



Niemela et al. |?fl(| have performed an extensive study 
of quasi- ID-rings with 2 to 4 interacting electrons. They 
used a homogeneous magnetic field and neglected the 
Zeeman splitting. The results for two and three inter- 
acting electrons show periodicities $o/2 and $o/3, re- 
spectively, i.e. consistent with the $>o/N periodicity for 
a narrow ring. In the case of four electrons Niemela et 
al. studied in addition to the 1/r Coulomb interaction 
also a (5-function interaction (VoS(r)) in the case of an in- 
finitely narrow ring. It is then interesting to compare the 
results of these two models, as shown in Figure El F° r 
the ^-function interaction the results are expected to be 
the same as for the Hubbard model with an infinite num- 
ber of lattice sites. Indeed, the results show the change 
of periodicity from $ first to $o/2 and then to $o/4, 
when the strength of the (5-function interaction (Vo) is 
increased. The results of Fig. El compare well with the 
results of the Hubbard model with only 8 sites, i.e. those 
shown in Fig. El for U = 0, U = 15 and U = 50. 
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The right hand side of Fig. ^] shows the energy levels 
for a quasi-lD-ring with 1/r- interaction. The spectra are 
essentially the same, the only difference being a slight 
upward shift of the energy when the flux changes from 
to 1. This is due to the harmonic repulsion of the flux- 
dependent effective potential, Eq. I|H9|I. caused by the 
homogeneous magnetic field. Note the similarities of the 
spectra shown in Figs. Elk and c. 



XII. PERIODICITY OF PERSISTENT 
CURRENT 

A. Strictly ID rings: Spectrum of rigid rotation 

The many-electron excitation spectrum for electrons 
interacting with the infinitely strong 5-function interac- 
tion was studied using the Hubbard model in Sec lIXI see 
Fig. [3 As mentioned, this spectrum can be constructed 
from the Bcthe ansatz. In the low-energy part of the 
spectrum (close to the yrast line) the levels consist only 
of the compact states of the quantum numbers Ij, i.e. 
there is no vibrational energy. It is instructive to use 
the continuum limit (L/N — > oo)) of the Bethe ansatz 
solution for the t-model to compare these (yrast) energy 
levels as a function of the flux $ to those of single electron 
states in a ring of radius R. One finds for a many-electron 
state 

E ^ $) = 2^R 2 ( M - hjk) 2 

(70) 

and for a single electron state 

h 2 ( $ \ 2 
< m ^=2^[ m -hTe) (71) 

where we have now written the flux quantum as $o = 
h/e. We notice that the many-electron states are iden- 
tical to the single electron levels with the electron mass 
and charge, m e and e, replaced by the total mass and to- 
tal charge of all the electrons, Nm e and Ne. Indeed the 
strongly correlated electron system with infinitely strong 
delta function interaction behaves as a rigidly rotating 
single particle. We should note that in the spinless case 
the number p of Eq. (|39[) is always zero and consequently 
one recovers the noninteracting case (for spinless elec- 
trons the (5-function interaction does not have any mean- 
ing due to the Pauli exclusion principle). The rigid rota- 
tion of the electron system then leads always to the $o 
periodicity as discussed already in Sec lIIII 

We have learned in Sees. [X] and IXII on the basis of 
numerical solutions for narrow quasi-lD rings, that elec- 
trons with normal 1/r interactions also produce simi- 
lar rotational spectrum. Moreover, the solution of the 
Calogero-Sutherland model shows rigorously that a simi- 
lar spectrum is observed for l/r 2 -interaction. Both these 
interactions have the property that they are infinitely 



strong at the contact, preventing the electron to pass each 
other. It seems obvious that any repulsive interaction 
which is infinitely strong (in such a way that electrons 
with opposite spin are not allowed at the same point) 
produces for the strictly ID ring the same spectrum of 
rigid rotation. 

The yrast spectrum is qualitatively similar for all elec- 
tron numbers. In particular, in experimentally deter- 
mined spectra, the number of electrons in the ring can 
be seen as a qualitative change of the spectrum only by 
observing the vibrational bands. (The yrast energy alone 
determines the number of electrons in a narrow ring only 
if the flux is quantitatively determined.) A way of esti- 
mating the number of particles is to count the number 
of (purely rotational) states below the onset of the first 
vibrational band, as can be seen from the following ar- 
gument: By considering noninteracting spinless electrons 
(which have exactly the same energy levels, though not 
all of them, as spinful electrons) we know from Eq.JSJ) 
that the first vibrational level is the first noncompact 
state, having an excitation energy of about Nti 2 /2m e R 2 . 
This equals the E(M = N, $ = 0) energy of the yrast 
band, Ea. (|70|) . For N particles there are thus about N 
purely rotational states below the first vibrational band. 
Figure [7| demonstrates that this is true for N = 4. 



B. Periodicity change in quasi- ID rings 

In this section we will concentrate on the lowest en- 
ergy state (the yrast state) and study in more detail its 
periodicity in quasi-lD rings, where the electrons are al- 
lowed to pass each other. (The periodicity of the persis- 
tent current at zero temperature is the same as that of 
the ground state energy). As discussed in Sec. IIX Fl the 
(strictly ID) Hubbard model suggests that (i) for spinless 
electrons the periodicity is always the flux quantum $o 
and (ii) for electrons with spin the periodicity changes 
from $o first to &o/2 and then to &o/N when the inter- 
action U increases from zero to infinity, as illustrated in 
Fig. [5] This change of periodicity in the Hubbard model 
was first studied by Yu and Fowler [55j and Kusmartsev 
et al.fr^. 

The change in periodicity can be traced back to the 
notion that the Hubbard model with empty sites can be 
quite accurately described by the t J c ff -model with an ef- 
fective exchange coupling J e g which depends not only 
on U but also on the number of empty sites, as demon- 
strated in Figs. ITrjlandlTTl When J c g is small, the $o/-^ 
periodicity is observed, while for large J e ff the $0 pe- 
riodicity is found (for even numbers of electrons). The 
periodicity of $o/2 results from the fact that the solu- 
tion of the Heisenberg Hamiltonian has two close lying 
states corresponding to angular momenta M — N/2 and 
M = 0, while the states corresponding to other angu- 
lar momenta are clearly higher in energy. The situation 
is demonstrated in Fig. 1161 where we show the energy 
spectrum for the Heisenberg model for an eight electron 
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FIG. 16: Schematic illustration of the &o/N and $o/2 pe- 
riodicities of a ring with eight electrons. The black points 
show the energy spectra with two different values of J. The 
parabolas drawn at each point indicate the change of the ro- 
tational energy level as a function of the magnetic flux. The 
lowest envelope of the overlapping parabolas gives the ground 
state energy as a function of the flux. 



ring. The energy spectrum for the total Hamiltonian as a 
function of the flux can be estimated simply by drawing 
a parabola at each energy level as shown in the figure. 
Now depending on the energy splitting of the Heisenberg 
model, i.e. on J e ff, the resulting lowest energy state as 
a function of flux can have periodicities, $>o/N, < &o/2, or 

Deo et ai.ji^l have studied the periodicity change using 
the model Hamiltonian Eq. l|t)8|) which is known to give 
good agreement with exact diagonalization results of the 
quasi- ID continuum Hamiltonian. This model Hamilto- 
nian has the same flux dependence as Ea. (|70fl . which was 
derived from the Hubbard model. It is thus natural that 
the same periodicity change is observed. In both cases 
the flux only affects the rotational state of the system by 
changing its energy as M 2 — > (M — N&/$ ) 2 , but does 
not change the internal structure of the many-body state 
in question, i.e. all interparticle correlations remain the 
same. 

We can look at the periodicity change in more detail 
by making a Fourier analysis of the lowest energy state. 
We have done this using the model Hamiltonian (|68H . 
A parameter determining the relative weights of the en- 
ergy states of the two separable parts of the Hamilto- 
nian is then JJ, i.e. the product of the moment of iner- 
tia (J ~ Nm e R ) and the Heisenberg coupling parame- 
ter (in the Hubbard model the corresponding parameter 
would be J e ff). Figure lTTI shows the three most important 
Fourier components as a function of J I for rings with 6, 
7 and 8 electrons. For 6 and 8 electrons we see clearly 
the period changes. In the case of an odd number of 
electrons, the solution of the Heisenberg model is quali- 
tatively different: There are two degenerate minima cor- 
responding to different angular momenta. Consequently, 
the period stays always more important than the 



FIG. 17: Three largest Fourier components of the periodic 
ground state energy. J I is the product of the effective Heisen- 
berg coupling and the moment of inertia. The dotted, dashed 
and solid lines are the Fourier components corresponding to 
the periodicities <E>o, $0/2 and <&o/N, respectively. 

period $ |I3- 

We have learned that the periodicity change as a func- 
tion of the effective width of the quasi-one-dimensional 
continuum ring is similar to that in a strictly ID Hubbard 
ring as a function of U, suggesting that the finite U in 
the Hubbard model mimics the finite width of a contin- 
uum ring. This similarity can be understood as follows. 
With infinite U the electrons are forbidden to pass each 
other. This situation is similar to an infinitely narrow 
continuum ring with 1/r interaction between the elec- 
trons. When U gets smaller, the electrons (with opposite 
spin) can hop over each other, the better the smaller U 
is. Naturally, in the continuum ring the electrons are al- 
lowed to pass each other if the ring has a finite width. 
Decreasing U in the Hubbard model thus corresponds to 
making the continuum ring wider. 

Note that in order to see the periodicity in the Hubbard 
model, one has to consider systems where the number of 
electrons TV is smaller than the number of lattice sites L, 
i.e. empty sites are needed for 'free rotation' of the ring. 
The resulting periodicity depends both on the number of 
empty sites and on the on-site energy U. 



XIII. OTHER MANY-BODY APPROACHES 
FOR QUASI-1D CONTINUUM RINGS 

Quantum monte carlo 

A class of powerful tools to study the low energy states 
of many- body quantum systems are based on the Monte 
Carlo method|8l]]. These methods have been extensively 
applied also for quantum dots[H HI H HI HI and 
quantum rings |87U8S| with a few electrons. Monte Carlo 
approaches are most suitable for studying either the 
ground state properties (variational Monte Carlo and dif- 
fusion Monte Carlo) or average finite temperature prop- 
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erties (path integral Monte Carlo), and have thus not 
been able to produce the detailed spectral properties with 
the accuracy of the CI method. It also seems that even 
finding the correct ground state is not straightforward by 
using Quantum Monte Carlo 89| 

Pederiva et al. 88] have used so-called fixed-node dif- 
fusion Monte Carlo for studying six electron quantum 
rings. They calculated also the lowest excited states for 
M = 0. They found the ground state to have 5 = and 
the first excited state 5 — 2, in agreement with the CI 
calculations and the Hubbard model, while for the sec- 
ond and third excited states they obtained the 5 = 1 and 
5 = states in opposite order as compared to the CI cal- 
culations. Nevertheless, as the authors note, the energy 
differences in the narrow rings are so small that their dif- 
ference starts to be within the statistical accuracy of the 
Monte Carlo method. 

The Monte Carlo studies for small quantum dots and 
rings show that while these methods can predict accu- 
rately the ground state energy, they are not yet capable 
to give reliably the salient features of the many-body 
spectrum. 



Local density approximation 

The density functional Kohn-Sham method is another 
approach mainy suitable for the determination of the 
ground state structure. In applications to quantum dots 
and rings (for a review see |37|). the local spin-density 
approximation (LSDA) is usually made and the system 
is assumed to be strictly two-dimensional. Generally, the 
Kohn-Sham method is a 'mean field' method, where the 
electron-electron correlation is hidden in an effective sin- 
gle particle potential. This causes the interesting feature 
that the mean field can exhibit symmetry breaking and 
the total electron and spin densities can reveal the inter- 
nal symmetry of the ground state, for example the inter- 
nal shape of nucleus or atomic clus ter|90l or the static 
spin-density wave in a quantum dot|9lj. Indeed, in ap- 
plications to quantum rings, the LSDA indicated the lo- 
calization of electrons in an antiferromagnetic ring|7 1|. 
Nevertheless, we should add that although the LSDA 
can often eludicate the internal structure of a rotating 
system, the method is not foolproof: In some cases, for 
example in rings or dots with high enough electron den- 
sity, it will not break the symmetry. 

Systematic studies of quantum rings in terms of density 
functional methods have been performed in several pa- 
pers (92 , 93, 94j,|95j, and comparisons with 'exact' many- 
body methods show that the LSDA gives accurately the 
ground state energy 88] . The LSDA has been extended 
to so-called current-spin-density functional theory (CS - 
DFT) which can take into account the gauge field |96|. 
Viefers et al. have applied the CSDFT for studying the 
persistent current, i.e. the ground state energy as a func- 
tion of the magnetic field, in small quasi-lD quantum 
rings. For a four electron ring (with r s — 2.5, Cf = 10) 



they found the yrast line consisting of two states: At zero 
flux (and at <E> = <&o) the ground state had 5=1 while 
around $ = $0/2 the ground state had 5 = 0. These 
results are in agreement with the CI calculations and the 
results of the Hubbard model for four electrons. Similar 
agreement was found for a six electron ring. The spin- 
densities showed a clear localization of electrons in an 
antiferromagnetic ring. 

Density functional theory has the same problem as 
quantum Monte Carlo in that the determination of 
excited states is not straightforward (although time- 
dependent current-spin-density-functional th eory can 
provide some information on excitations[92, lU 
Consequently, it is not possible to construct the com- 
plete excitation spectrum as by using the brute force CI 
method. 



XIV. RELATION TO LUTTINGER LIQUID 

Infinitely long one-dimensional systems are often stud- 
ied as Luttinger liquids [9Hl9S| (for reviews see 
The speciality of the strictly ID systems arises from the 
fact that the Fermi surface consists of only two points 
(±&f). This leads to a Peierls instability[99| and a break- 
down of the Fermi liquid theory in a strictly ID system. 
Important low energy excitations will then be collective, 
of bosonic nature, and have a linear dispersion relation. 
The Luttinger liquid also exhibits so-called charge-spin 
separation: The spin and charge excitations can move 
with different velocities. In addition to studying the low 
energy excitations, the Luttinger model has been exten- 
sively used for studying correlation functions 0,0- 

It has been shown that the tJ-model is a Luttinger 
liquid 6[. In the limit of a narrow ring with many elec- 
trons the spectral properties of the quantum rings must 
then approach those of a Luttinger liquid. We will now 
demonstrate that the many-body spectra of quantum 
rings are consistent with the properties of the Luttinger 
liquid. We do this only by qualitative considerations. In 
an infinitely long ID system, the low energy single par- 
ticle excitations (of free fermions) are restricted to have 
a momentum change of q ~ due to the fact that the 
Fermi surface is a point. In the Luttinger model it is pre- 
cisely these excitations that lead to the bosonization and 
collective plasmon excitations 6] with a linear dispersion 
relation. In the case of a finite ring these single parti- 
cle excitations are just those described in Fig. ^ where 
one of the last electrons is excited from the compact state 
(or similarly in the Bethe ansatz solution of the i-model) . 
Now, it is exactly these single particle excitations which 
lead to the vibrational model, i.e. longitudinal acoustic 
phonons in the limit of a long ring (see Fig. 0J, which 
have a linear dispersion relation (for small q). The ex- 
citation spectrum is in qualitative agreement with the 
prediction of the Luttinger model already in the smallest 
rings. 

Casting the Hamiltonian explicitly into charge de- 
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pendent parts (rotations and vibrations) and a spin- i 
dependent part (Heisenberg Hamiltonian), as in Eas. (|13|) 
and (fr!r)f is equivalent to the charge-spin separation in the 08 
Luttinger model. In an infinite system this can be explic- 
itly done for the half- full Hubbard model (L — N). We 
have demonstrated in Sec. IIX Fl that the spin degrees of | 
freedom can be described with a good accuracy with the s 
Heisenberg model in a much larger variety of quasi-lD 2. 04 

rings. 



XV. PAIR CORRELATION 

The internal structure of a many-body electron state, 
especially the possible localization to a Wigner molecule, 
can be studied by examining the correlation functions. 
We have done this already for non-interacting electrons 
in Fig. |3 where we used the iV-particle correlation func- 
tion for identifying the vibrational states. The iV-particle 
correlation function is just the square of the normalized 
many-body wave function. A related analysis of the max- 
imum TV-particle correlation for the Hubbard model, in 
Table HI also revealed the internal stucture of the vibra- 
tional states. 

The pair correlation function is frequently used to 
study the internal structure of a many-body state. In 
one-dimensional systems the pair-correlations have the 
property that they decay with distance as l/r a j§| . Con- 
sequently, the ID electron system does not have a Wigner 
crystal with true long-range order. The same is true for 
the spin-density oscillations: For example in the antifer- 
romagnetic Heisenberg model in ID the spin-spin corre- 
lation decays as 1/r, r being the distance between the 
electrons. 

In finite ID rings, with only a few electrons, the pair 
correlation function is even less informative. The rea- 
son is again the Pauli exclusion principle, which prevents 
electrons with the same spin to be at the same site. This 
means that within a short distance, the pair correlation 
functions are quite insensitive on the electron-electron in- 
teraction. This is demonstrated in Fig. 1 181 where we show 
the calculated p air c orrelation function for two different 
continuum rings [Toof and compare them to the pair cor- 
relation of the Heisenberg model. We have also studied 
the pair correlation for four electrons using the i-model 
and found that the correlation is independent of the num- 
ber of empty sites, as expected from the notion (Section 
IIX F|) that the spin-spin correlation is determined from 
the Heisenberg model, whatever the number of empty 
sites. 

The effect of temperature on the pair correlation func- 
tion has been studied by Borrmann and Harting|87| us- 
ing quantum Monte Carlo. Koskinen et aL^?J used the 
model Hamiltonian (|66f) to determine the temperature 
dependence by calculating separately the pair correla- 
tion for each quantum state. Both methods agree in the 
fact that the correlations between the electrons vanish as 
soon as the temperature exceeds the first excited state of 
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FIG. 18: Pair correlation function for six electrons in a 
quasi-one-dimensional continuum ring, calculated with the CI 
method. Solid and long-dashed lines: |J.-correlation and TT _ 
correlation, respectively, for a narrow ring with r s — 2 and 
C'f = 25. Dashed and dotted lines: ^-correlation and It- 
correlation, respectively, for a wider ring with r s = 2 and 
Cf = 4. The filled and open circles show the correlations for 
a six electron antiferromagnetic Heisenberg ring. The func- 
tions are normalized so that the integral of the ||-correlation 
is 3 and that of the j j-correlation is 2. 

the system. 

XVI. INTERACTION OF THE SPIN WITH THE 
MAGNETIC FIELD: THE ZEEMAN EFFECT 

Throughout most of this paper we have assumed that 
the magnetic flux is confined inside the ring so that 
the electrons move in a field-free region. Experimen- 
tally, however, it might be difficult to produce a situation 
where the magnetic field is zero at the ring site (or the 
effective Lande factor is zero). It is then important to 
consider also the Zeeman effect when comparing theory 
with experiments. The interaction between the electron 
spin and the magnetic field adds to the Hamiltonian a 
term fiB9oS z B. In the case of a narrow ring it is benefi- 
cial to write this as 01 

ff^^ r )gS z , (72) 

where we have written the field with help of the flux, ring 
radius, and an effective Lande factor g. For example, in 
the case of a homogeneous magnetic field and an ideal ID 
ring $ = irR 2 B, and g = go- The advantage of writing 
the Zeeman part of the Hamiltonian as above is that we 
can study continuously the change in the spectra when 
we move from the case 11140 . no Zeeman effect, to a homo- 
geneous magnetic field simply by changing g. Moreover, 
the fact that the effective Lande factor in semicondictors 
differs from that of the free electron can also be taken 
into account by just changing g. 

The periodicity of the persistent current, or lowest en- 
ergy state, as a function of flux is caused by the ground 
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FIG. 19: Angular momentum of the ground state of a six 
electron ring as a function of the flux and the effective Lande 
factor g. The angular momentum is shown as numbers and it 
increases with the darkness of the grayscale. The uppermost 
panel is for J I = 6 in the model Hamiltonian l|66fl , the center 
panel for JI — 1.8 and the lowest panel for JI — 0. 



state energy jumping from one angular momentum state 
to another. If all possible angular momenta are visited 
(Jeff is small in the models discussed in Sec. IIX F|) . the 
periodicity is &o/N; if angular momenta M = 0, N/2, N 
etc. (or N/2, 3N/2 etc.) are visited the periodicity is 
<i>o/2, and if only angular momenta M — 0, TV, 2iV etc. 
are visited the periodicity is <I>o- I n order to determine 
the overall periodicity it is thus sufficient to examine the 
angular momentum values of the lowest energy state as 
a function of the flux. 

Figure IT01 shows the angular momentum of the ground 
state as a function of the flux and the effective Lande fac- 
tor discussed above. The results are shown for three dif- 
ferent values of the parameter J e ff. When the Lande fac- 
tor approaches the free electron value g — 2 and the flux 
is large, the periodicity always becomes $o- The reason is 



the large Zeeman effect which makes the maximum spin 
state S z — N/2 the lowest energy state. This state has 
always the periodicity $0 as shown in Fig. In a wide 
ring (a) where the periodicity is $0 for g = it changes 
first to 4>o/2 and then again to <&o when g increases. This 
is due to the fact that (in a ring of six electrons) angu- 
lar momenta M = 0, 6, • • • are the ground states for the 
nonpolarized case while for the polarized case the ground 
state angular momenta are M = 3, 9, • • • . In the transi- 
tion region a periodicity <&o/2 is observed. In the narrow 
ring (c) the periodicity changes smoothly from &o/N to 
$0 when g is increased from zero. 



XVII. EFFECT OF AN IMPURITY 

There has been extensive research on the effects of im- 
purities on the p ersistent current. Already in 1988, Che- 
ung et al. |l0l| used a simple tight binding model for 
studying the effect of disorder in ID rings and found 
the decrease of the persistent curren t wit h increasing 
disorder. Chakraborty and Pietilainen |l02j used the CI 
method for studying the effect of an impurity in a po- 
larized electron ring. The basic results of the included 
Gaussian impurity potential were to lift some of the de- 
generacies of the energy levels and decrease the persistent 
cu rren t. Similar findings were observed by Halonen et 
al. 103] in the case where the the spin was included. The 
suppression of the persistent current due to a gaussian 
impurity was also computed by Viefers et aL[24( using 
density functional theory, thus including the effects of 
both spin, realistic i nter actions and finite width of the 
ring. Eckle et a/.j&iJES] have used the Hubbard model 
to study the effect of coupling the quantum ring to a 
quantum dot at Kondo resonance. They have shown that 
at certain values of the flux the problem can be solved 
exactly using the Bethe ansatz. 

The effect of an impurity is easy to study in the Hub- 
bard model, where in the simplest case we can just add a 
repulsive or attractive potential in one of the lattice sites. 
Figure [3U| shows how the flux dependence of the energy 
levels changes when a repulsive external potential is in- 
troduced in one of the lattice sites. The effect is to open 
gaps in the spectrum. The impurity potential splits the 
degeneracy of angular momenta (modulo 4) and makes 
all the energy levels oscillate with a period of <I>o instead 
of the period of NL&o seen in the impurity-free case. 
Note, however, that the levels corresponding to different 
angular momenta still cross and the lowest energy state 
has the same period as without the impurity. 

In the above example, the energy of the lattice site, 
say site number 1, was changed by introducing a term 
A(nif + nil) m the normal Hubbard Hamiltonian, Eq. 
(J22J. Figure |2*T1 shows how the lowest energy state as 
function of flux changes when A is varied from -20 to 20 in 
steps of 2. For a relatively small value of U (=20) and an 
attractive impurity, the variation of the ground state en- 
ergy as a function of the flux is nearly independent of the 
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FIG. 20: Effect of an impurity on the energy levels of a six- 
site Hubbard ring with four electrons. The left panel shows 
the energy levels as a function of the flux (in units of $o) for 
an impurity free ring. The middle panel shows the evolution 
of the energy levels at zero flux as a function of an impurity 
potential A at one of the lattice sites. The right panel shows 
the energy levels as a function of flux in the case where A = 1. 
U — 1000 in all cases. 
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FIG. 21: Effect of an impurity on the lowest energy state in 
a Hubbard ring with L — 6 and N — 4. The left panel is for 
U = 20 and the right panel for U = 200. The thick line is the 
result without the impurity and the other lines correspond 
to rings where the energy of one lattice site is increased (see 
text) by 2, 4, 6, etc. (dotted lines above the thick line) or 
decreased by 2, 4, 6 etc. (dotted lines below the thick line). 
All lines have been shifted vertically so that they appear to 
be equally spaced. The energy scale in each line is the same, 
but it is twice as large for U = 200 as compared to U = 20. 
The flux is in units of <E>o = h/e. 
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FIG. 22: Energy spectra as a function of flux for four elec- 
tron clusters. Upper panel: Hubbard ring with five sites and 
finite (7 = 7. Middle panel: "Quasi-one-dimensional" f-model 
ring. Lowest panel: Quasi-one-dimensional t-model ring with 
a defect. The geometries in each case are shown at the left. 



impurity potential, while a positive (repulsive) impurity 
potential reduces the amplitude of t he oscilla tion, as ex- 
pected from all earlier studies[7(10[iol[iol]. Since the 
ground state persistent current is essentially the deriva- 
tive of the energy wrt. flux, see Ea. (|24(l . this directly 
reflects the suppression of the current. In the case of 
a relatively large U (in the present example U — 200), 
we see that the impurity potential changes the original 
periodicity <J>o/iV to $o/2. 

In the case of small systems the Hubbard model can 
easily be to used in studying also other kinds of impu- 
rities. For example, two impurities or even random po- 
tentials at lattice sites have qualitatively the same ef- 
fect on a Hubbard ring as a single impurity, as shown in 
Fig. 1201 The Hubbard model may furthermore be used 
to study the coupling of the ring to a quantum dot[Hflj. 
as mentioned earlier. In Figure we present another 
example demonstrating how the Hubbard model can be 
extended to quasi-lD rings with impurities. In the upper- 
most panel of the figure we show the energy spectrum of 
a ID Hubbard ring with four electrons and a finite U = 7. 
The middle panel represents a quasi-YD ring with infinite 
U (t- model). Comparison of these two models shows that 
a finite U in a strictly ID Hubbard model actually mim- 
ics well the quasi-one-dimensionality of a more realistic 
ring, as noticed already in comparing the Hubbard model 
with the continuum models. The lowest panel of Fig. l2*2l 
shows an 'impurity' or a narrow neck in the ring. Its 
effect is to open gaps in the excitation spectrum. Qual- 
itatively, the effect is the same as observed above for a 
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ID Hubbard ring with an impurity potential. 

In many experiments, the persistent current is mea- 
sured for a collection of rings, possibly consisting of sev- 
eral propagating channels (radial modes). The persis- 
tent currents in such rings are not determined by quan- 
tum mechanical eigenenergies of a single ring, but are 
affected by sever al complication s like disorder and en- 
semble averaging |l05l Il06. 107]. The electron-electron 
interactions s eem to play a crucial role also in the dis- 
ordered rings |l08| . Miiller-Groeling and Weidenmiillcr 
have shown that interactions effectively counte ract the 
impurity suppression of the persistent current fl09| . in 
qualitative agreement with our Hubbard model results. 

XVIII. SUMMARY 

In this paper we have attempted to present a compre- 
hensive review of the physics of few-electron quantum 
rings, with particular focus on their energy spectra and 
the periodicity of the persistent current. We compared 
various analytical and numerical theoretical approaches 
which fall into two main classes - lattice models on the 
one hand and continuum models on the other - and tried 
to clarify the connections between them. The main mes- 
sage is that all the different approaches give essentially 
the same results for the spectra and persistent currents 



(provided one takes the continuum limit of the discrete 
models to compare them to the continuum ones). The 
essential physics is captured by a simple model Hamil- 
tonian which describes the many-body energy as a com- 
bination of rigid rotation and internal vibrations of the 
ring 'molecule', plus a Heisenberg term which determines 
the spin dynamics. In the spectra, the vibrational excita- 
tions can be seen as higher bands, while the lowest (yrast) 
band is purely rotational. Its periodicity with respect to 
the flux is always $0 m the case of spinless (polarized) 
electrons, while for a clean ring of nonpolarized electrons, 
the periodicity changes from <S> via <&o/2 to &o/N as the 
ring get narrower or, in the language of the Hubbard 
model, as the interaction strength U increases (the ef- 
fectvie Heisenberg coupling J decreases). Impurities in 
the ring may change the periodicity (to $0) even in the 
nonpolarized case, and moreover suppress the persistent 
current. 
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